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Abstract 

In terms of multiple scattering picture, we calculate the double scattering effect in the transverse 
momentum distribution of J/tp photoproduction. Applying the generalized factorization theorem, we 
find that the contributions from double scattering can be expressed in terms of twist-4 nuclear parton 
correlation functions, which is the same as that used to explain the nuclear dependence in di-jet mo- 
mentum imbalance and in direct photon production. Using the known information on the twist-4 parton 
correlation functions, we estimate that the double scattering contributes a small suppression in J/ip pho- 
toproduction. In the analysis we only take into account the leading order in the small velocity expansion 
for the nonperturbative parts related to the quarkonium. 
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1 Introduction 



Heavy quarkonium production in high-energy collisions affords us a variety of insights into the underlying 
dynamics of the strong interaction. Especially in nucleus- nucleus collisions, we may find that a suppression 
of ,// -0 production can serve as a clear signature of the existence of the deconfining phase of QCD [Q . This 
suppression effect was observed by NA38 collaboration later ||. However, successive observation || pointed 
out that in proton-nucleus collisions, in which there is no QGP formed, such suppression could also exist. 
There have been a number of attempts to explain it, such as pre-resonance absorption, gluon shadowing, 
hadronic co-mover absorption and the energy loss model ^, [|, |[ @, ||, [§, etc. To understand the J/ip 
suppression mechanism clearly, its production by scattering in different nuclear facilities must be studied 
carefully. 

In principle, the J/ip state can be described in a Fock state decomposition 

| J/0 >= 0(l)|cc( 3 sf 3 ) > +0(v)\cc?pf ) )g >+-.., (1) 

where 2S + l i/j^ characterizes the quantum state of the cc pair with color-singlet (CI) or color-octet (C8) 
respectively. This expression is valid for the non-relativistic QCD(NRQCD) framework and the coefficients 
of each component depend on v, the small relative velocity between quark and antiquark. The charmonium 
production can be divided into two steps. The first step is the production of cc pair. The cc pair can be 
either (cc)i or (cc)s, and are produced perturbatively and almost instantaneously, with a formation time 
Tf sa (2m c ) _1 = 0.07 fm in the cc rest frame. The second step is the formation of a physical state of J/ip from 
cc pair, that need longer time. The second step is a nonperturbative process, and in different environments 
such as vacuum, hadron matter and QGP, the original produced cc pairs undergo different transition rules 
to final observable J/ip. All these processes have been studied carefully before. However, besides this 
nonperturbative effect related to the transitions there is another nonperturbative effect if there are hadrons 
or nucleus in the initial states. This nonperturbative effect can be analyzed with a twist expansion method 
PI - Ref. (l^] gives the effect of the order of next-to- leading twist in inclusive photoproduction of 
quarkonium. In this paper we will use factorization at higher twist to describe multiple scattering effect for 
J/ip production in photo- nucleus collision, and our result can be extended to study the corresponding cases 
of hardon-nucleus and nucleus-nucleus collisions. 

In terms of factorization at higher twist, Luo, Qiu and Sterman (LQS) |l3L [l4| have developed a con- 
sistent treatment of multiple scattering at partonic level. They showed that nuclear effect may be brought 
directly into the scattering formalism of QCD, by treating it as a factorizable nonleading-order correction 
to hard scattering, and nuclear enhancement appears in this context as a property of multi-parton matrix 
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elements. In |13| they derive the anomalous nuclear dependence of jet cross sections in elastic scattering 
and photoproduction in terms of twist-four parton distributions in nuclei. In Jl5| they expressed the nuclear 
dependence of di-jet moment imbalance in photo-nucleus collisions and estimated the size of the relevant 
twist-4 parton distributions by using the Fermi Lab E683 data. Recently Guo applied this method to 
Drell-Yan process, and calculated its transverse distribution in hadron-nucleus collisions. 

All above processes are totally inclusive, and the extension to our case is not quite straightforward. The 
process we consider is seminclusive, where a cc pair is observed indirectly. We neglect the higher order v 2 

contributions and consider that the formation of J/ip is totally evoluted from the leading order cc( 3 s[^) 
pair. This formation starts with a CI or C8 cc pair produced from initial hard collision between photon and 
a parton from nucleus, and then on its way out of the nucleus, it undergoes a soft scattering with another 

parton from the nucleus to form a color single cc^S^) pair. The soft scattering for cc pair can be done 
by either its quark or its antiquark picking up the soft gluon from nucleus. Our paper is divided into four 
parts. In section II, we derived the general formalism and the result of single scattering. In section III, the 
double scattering contributions are derived. In the last section, we present the numerical results and give a 
brief discussion. 



2 General formalism and the single scattering result 

To study the nuclear effect in transverse momentum distribution, we consider the differential cross section 
for the inclusive process: 

70/) + A(p) J/V(Q) + X. (2) 

The momenta of particles are given in brackets, p' is the momentum of incoming photon, p is the momentum 
per nucleon for nuclear target, and Q is the momentum of outgoing J/^p. In terms of contributions from 
multiple scattering, we expand the differential cross section as 



dQ-^dy dQj_dy dQ]_dy 

The superscript "S" and "D" denote contributions from the "single scattering" and "double scattering", 
respectively, and "• • •" represents contributions from even higher multiple scattering. We choose a frame in 
which the nucleus A moves in the ^-direction and the photon in the opposite direction, Q± is the transverse 
momentum of the J/ip, and y is its rapidity 

y = 1 ln Qa + Qz (4) 
2 Qq — Q z 
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Figure 1: A Graphical representation for single scattering process in photon-nucleus collision. 



In our approximation the quarkonium mass is twice of quark mass m c . We define quantities s = (p+p 1 ) 2 , 
t = (p — Q) 2 and u = (p' — Q) 2 . For performing the analysis it is convenient to use light-cone coordinate 
system. In this system the photon carries the momentum p ,f * — (0,p'~ , OjJ, and p^ = (p + ,0, OjJ. We 
introduce in this frame two vectors n and h and a tensor d^: 



= (0,1, 0_l), = (1,0, 0_l), 



d^ = cf v - n^h v ~ n^ft". (5) 



We will work in the light-cone gauge n ■ A(x) = A + (x) = 0, where A fl (x) = A a ' tl {x)T a is the gluon field. 



In order to compare with experimental data, we calculate the ratio of total differential cross section and 
single scattering contribution, 



= da jA / dgS A ^ da° A I da^ A 
~ dQ\dy I dQ\dy~ dQ\dy / dQ\dy 



To analyze the effect of single scattering we need to consider the diagram in Fig.l, in which the up box 
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represents the nonperturbative part related to the nucleus A, the down box contains the perturbative part 
for production of the cc pair at the leading order of coupling constant and the transition of the pair into J/ij), 
they are illustrated in the Feynman diagram given in Fig. 2. The contribution from Fig.l can be written as: 

2sa " A = J !0fS;\%{k) J d A ze ik - z <p A \A a ^(z)A a ^(0)\ PA >, (7) 

where S?*?£(k) corresponds to the lower part in Fig.l, its leading order part is shown in Fig. 2, where the 
left and right box include all possible tree Feynman diagrams with the external partons shown in the figures. 
In the center of mass frame of high energy collision, all partons inside the nucleus are parallel to each other, 
along the momentum of the nucleus. In order to get the leading order contribution, we make the collinear 
expansion, 



Afc" 

k—xp 



2 dk a dkP 



Afc Q Afc /3 + • • • , (8) 

k—xp 



Where Afc = k — xp. We only keep the first term S"lul (xp), which is independent of k and k±, so we can 
perform the integration over k~ and k± directly, 

2*4, = J ^S™%{xp)^ J dz~e^ +z - < PA \A a ^(z)A a ^(0)\ PA > 

where we have used the matrix element relation JL2| 

< PA \A a ^(z)A a ^(0)\ PA >= Ld^S^ <p A \A a± (z)A a ± (0)\p A > . (10) 
The function f g / A (x) is defined as 

fg/A( x ) = / dz-e****'- < p A \A a± (z)A a A0)\p A > . (11) 

We consider the formation of J/ip from a cc pair at leading order in v 2 . According to the NRQCD 
factorization formalism [fL7| , the amplitude of J/tp production from partons can be factorized into 
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/J>2, 0,2 




Figure 2: The diagram corresponds to the lower part in Fig.l for J /if) production at lowest order partonic 
level. 



A(gj -> M = A( 91 -> cc( 3 s[ 1] )) pert . ® AiccCS^) . (12) 

nopert. 

Here we only keep the leading order contribution comes from color singlet cc( 3 S± ) pair. Contribution from 
color octet and other higher angular momentum cc pair states are suppressed at least v 2 , we neglect them 
here. cc^S^) pair production starts at 0{a 2 s ) via process 75 — > V J ffi whose Feynman diagram is shown in 
Fig. 3. Now S'SiSCxp) can be written as 

= \t (2) x L£ vX (k) ® RZ^ik) x X(V^ C ), (13) 

where i represents the average over the initial photon spin, 

r(2) = dQldy 1 u 

8?r s + t-m^ v s + t-m^ v ; 

represents the final two body phase space integration for J ftp and a gluon. Nonperturbative part M(j/j c ) in 
Eq.(|l3|) cannot be computed from first principles without resorting to lattice QCD, we neglect the low energy 
dependence of interactions, the nonperturbative amplitude only affect the overall normalization constant, 

M{4> c ) = — !- < 0\OlP( 3 Si)\0 >= constant. (15) 



In ( [l3| ) L°^ iuX and i2xt/*a represent the perturbative production amplitude the of cc( 3 <S'i 1 ' ) ) pair shown in 
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Figure 3: Diagrams represent the cc pair production amplitude including all partonic hard subprocess 
leading order in perturbative coupling constant. 



the left and right parts of Fig. 2 respectively. We have 



L l\u\( k ) = -^^ J u{qi]s 1 )L a ^ vX {k,q u q 2 )v{q2,. 



si,s 2 



V3 



^ux^q^q^PisAquQh) > (16) 



where L°^ iV ^{k 1 q%, q 2 ) represents the left box of Fig. 2 with external quark momentum qi and antiquark 
momentum q 2 , and its all tree level diagrams are illustrated in the diagrams of Fig. 3. The heavy quark and 
antiquark propagate nearly on -shell with a large combined Q = q\ + q 2 and a small relative momentum 
q = qi — q 2 . The spin quantum numbers for the cc pair are projected out by the sums over SU(2) Clebsch- 
Gordan coefficients < isi; ^s 2 \lS z >. The 4x4 matrix is 

As, = ^2 v(q 2 ; s 2 )u(qi; si) < ^si\^s 2 \SS z > 

sis 2 

_i ' -m c )i{Q,S z ){i 2 +m c ) (17) 



2 v / 2m c j 

up to 0(q 2 ) corrections. The project operators' dependence on q is irrelevant for S'-wave cc production. This 
indicates qi = q 2 = q, and Q = 2q, so 

P ls ^ = -Z}—(i- mc )^Q,S z ){i + m c ). (18) 
2v2m c 

The final perturbative result can be derived asfl2]| 

da s lA _ 327:aay Q m^\R(0)\ 2 f g/A (x) f P(§ - m\ f + t 2 (t 2 - m\ f + (s + t) 2 (s + t - m% f ' 



(19) 



dQ 2 ± dy 3s( S + i-m^) x [ (§ - m% ) 2 {t - m% ) 2 (§ + t) 2 

where 

s= (xp + p') 2 , t=(xp-Q) 2 , (20) 
and x is fixed according to the momentum conservation, 

x = -—^— T . (21) 
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Figure 4: Graphical representation of double scattering contribution in photon-nucleus collisions; a) 
diagram, b) and c) interference diagram. 
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3 The double scattering result 



Now we begin to consider the double scattering effect, at lowest order, the cross section corresponding to 
soft rescattering are shown in Fig. 4, its contribution can be written as 



2sda° A 



d 4 hd 4 k 2 d 4 k 3 

(27r) 4 (27r) 4 (27r) 4 "w/*/*/* 1 
< p A \A a ^{z 1 )A a ^{z 2 )A a3 ^(z 3 )A ai ^ i {Q)\p A > 



(22) 

The momenta, color- and Lorentz- indices are marked in Fig. 4. In Eq.(|2|) C®^^®*(ki, k 2 ,k 3 ) corresponds 
to hard production of cc pair and its transition to J/ip- To pick up the leading contributions to the nuclear 
enhancement, we expand the parton momenta at the following according their collinear values: 

fci = xip + fi, 

k 2 = x 2 p + n + k 2X , 

h = x 3 p + n + k 3J _. (23) 

The minus components of the fcj give even smaller contributions in C^^f^^ki, k 2 , k 3 ) and will be neglected. 
In addition, k%± dependence in C^^mI^ will not associated with A enhancement and we drop it as well. 
Later, we will show that x 2 and 2:3 can be fixed by poles as functions of k 2 ± (or k 3 ±) and Q±, and that they 
vanish when k 2X and k 3X go to zero. 

Once we have dropped the k\ X and k~ (i=l,2,3) dependence in C*^^^*, their integrals give S functions 
and allow us rewrite Eq.(p2) as 



2sda° A 



f p+dx! f p+dx 2 d 2 k 2± fp+dx 3 d 2 k 3± f 2 2 

J —J J —w? J 1 2 Z2± 3 Z3J 



C a ^llll\(k u k 2 ,k 3 )< PA \A^^(z 1 )A a ^ > . 

First of all, we expand the hard part C^^f^f 4 (ki, k 2 , k 3 ) in momenta about the collinear direction: 



(24) 



S~<ala 2 a 3 a4, 



dC*\%%%{x lP ,k 2l x 3 p 



dk p 2X 



K 2_L 



dC a ^Zl\(x lP ,x 2 p,k 3 ) 



kn , =0 



dk p 3X 



K 3_L 



k 3 ±=0! 



a2 /ya 1 a 2 a 3 a 4 , ( h, u \ 

° °MiM2M3M4 ( - a; lP' /C 2 J K3j 



k p k a 



(25) 



k2±=k 3 ±=Qi 



10 



In this discussion we will confine ourselves to experiments with unpolarized beams. Thus odd- twist contri- 
butions of the second and third terms in Eq.(^5|) vanish on taking a spin average. In addition, the leading 
term C^'^ l 3 3 a l f i (xip, x 2 p, x 3 p) does not contribution to A dependence because it is independent of transverse 

momenta. We did not list explicitly terms such as d 2 C/dk 2p and d 2 C/dk 2tT because, for reasons that will 
become clear below, they do not contribute physical double scattering. Therefore, double scattering effect 
comes entirely from the last single term of Eq. (p5[) 



2sda° A 



1 f (p^dxidxidxs f _ ( d 2 C%%%%\(x lP ,k 2 ,k 3 

2 J (2tt)3 J Zl H H { 



e ix 1 p + z 1 +ix 2 p + z 2 +ix 3 p + z 3 

d 2 k 2 ± d 2 k 3J _ 



(2tt) 2 (2tt) 2 



fc 2 \fc 3 \e lfc — +lfc — < p A \A^( Zl )A^(z 2 )A^(z 3 )A^(0)\p A > 



2 J (2tt)3 J 1 : - :;; I 



e ix 1 p + z 1 +ix 2 p + z 2 +ix 3 p + z a J ^ 



AaAZl) dz 2p dz 3rT Aa * (0) 



PA > 



(26) 



Because the values of x%,x 2 and x 3 will be fixed by 5 functions and poles of the hard part C*^^^° 4 4 , we 
move the exponentials e lXlP+z ± + lx 2p + z 2 +"3P + z 3 anc [ integrals inside the derivative of d 2 / dk 2 dk 3 , 



2sda° A 



1 /" (p + ) 3 dz 1 dz 2 dz 3 



(2tt)3 



< Pa 



< 4 (0) 



Pa > 



dx 1 dx 2 dx 3 e ix ^ + ^ +ix ^ + ^ +ix ^ + ^ 6«\%% a p \ (x lP , k 2 ,k 3 ) 



(27) 
(.28) 



fe2J_=fc3_L=0j_ 



According to the discussion in p2| [T3| , the dominant contribution of the four-gluon matrix element can be 
written as 

< PA \A%te)G%Pte)G£°(zt)A>£(p)\p A > 



1 V^V^ , m 



4x8x8 (p+) 2 



d^df j«i»45»2a 3 < ^ |^- L ( z -)G'; ± (z 2 -)G' b+± (z 3 -)yll(0)|p J 4 > (29) 



Therefore, the cross section for double scattering turns out to be 



D 1 1 f p + dz 1 dz 2 dz : 



(2tt) 3 



<p A \A^{z-)G b + ^2)G b+ ^{z 3 )Aimp A > F(z-,z-,z-), (30) 
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where 



(31) 

K2J_=«3X=0_L 

Equation (|3fJ) has been derived at lowest order in the hard scattering. As it stands, it is not gauge 
invariant. Gauge invariant is incorporated at higher orders and leading power by taking into account diagrams 
involving more fields n ■ A, which will form ordered exponentials between the physical fields shown in J30| ) , 
as discussed, for instance, in jul Such fields, however, can not correspond to physical rescattering, since 
they can be eliminated by a change of gauge. This is the reason why we neglected two derivatives with 
respect to ki and none with respect to kz , or vice versa. Such a term would involve only three physical fields 
instead of four. As we shall see below, four physical fields are required to produce nuclear effect. 

From Eq.@ and @ all integrals of £i)£2j#3 can now be done explicitly without knowing the de- 
tail of the multi-parton matrix elements. Now we are ready to discuss the calculation of the hard part 
C^^°|^. 4 4 (a;ip, k2, k^). As is shown in the diagram of Fig. 4, during the collision, firstly the photon strongly 
interacts with a gluon from one nucleon in nucleus. After a hard scattering, the outgoing cc pair interacts 
with a soft gluon from another nucleon in the nucleus to form a color-singlet cc( 3 Si) pair. So there are 
two partons from nucleus participate in the scattering process. The kinematics can only fix one parton 
momentum, thus we need to integrate over the other parton's momentum. In C^^^."^ {x\p, &2, ks) we will 
encounter one phase space S function and two virtual propagator poles involving x^s, which correspond to 
the zero momenta fraction partons. The leading contribution of the integration over the extra parton mo- 
mentum is given by the residues at these poles, which will fix xi and X3 to be functions of the transverse 
momentum Then in the limit kij_ — > the X2 and X3 will vanish. Therefore e lx 2P +z 2 and e lx 3P +z 3 W i\\ 
reduce to unity eventually. We may expect a significant contribution from the free integrals J dz^dz^. This 
will be the origin of the A dependence. 

Consider subprocess shown in Fig. 4, there are three four-momentum linking the partonic part and cor- 
responding four gluon matrix element. In inclusive DIS fll3| and Drell-Yan [ fl6| processes, the final results 
depend only on the real diagram with one extra gluon in each side of the cut, and the role of interference 
diagrams with both gluons in one side of the cut is to take care of the infrared sensitivities of the short 
distance hard parts. The process we are considering is seminclusive and the subsequent cc pair is in color- 
singlet state. This demands that for the real scattering subprocess shown in Fig. 4(a), the original cc pair 
produced from hard scattering of photon and one parton should be in color-octet state, while for the virtual 



y /)a 1 a 2 a 3 a i ( u iOr/' 2 r)' i3 ,7 A ' lA14 A f) 1 
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(a) (b) 

Figure 5: Graphical representation of the left amplitude L a ^J x in Eq.(|32|) for cc interacting with soft gluon, 
(a) shows the quark picks up the gluon, (b) shows the antiquark picks up the gluon 



subprocesses shown in Fig. 4(b) and Fig. 4(c), the corresponding cc pair should be in color-singlet state. Thus 
we must consider the two different types of physical rescattering separately, and as we shall see below, each 
of their contributions remains infrared finite in itself. This is quite different from the inclusive situations, in 
which divergences are canceled only after summing over all sorts of cut in a Feynman diagram. 

First of all, we consider contribution of real double scattering diagram illustrated in Fig. 4(a). By the 
similar approach which we have taken in deriving Eq.([l3|), we can write its hard part Cj^aMsX as 

Cl\Z% a ^ l \x lP , k 2 ,k 3 )p^p^ = ir( 2 > x L«l a J x (x lP , k 2 ,k 3 ) ® RZI'MP, k 2 ,h) x Mty c ), (32) 



Nonperturbative part A4(jJj c ) has been defined in Eq.(|l5[). The final state gluon quarkonium two particle 
phase space can be written as 

(2) dQ\dy 1 u - 2Q ■ k 2 ± + fcf ± . , , 

r = — Z T7 2 S(x 1 +x 2 + — r^)- 33 ) 

87r s + t- mi s + t — mfi, 



The represents of amplitudes L" 1 ^^ and R\ 3 J^ 4 are shown in Fig. 5 and Fig. 6 respectively. The cc pair 
can interacts with nucleon through one of its constituents quark or antiquark scattering with soft gluon. In 
case of the quark participating the interaction, as shown in Fig. 5(a), we have 

T a ia2 (a), , , x „ f ij- m c )j{Q, S z ){j + m c )j){i - x 2 j> - fe± + m c ) f , \ 
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f (4 - m c )({Q, S z ){i + m c )(2p ■ q - ± )T a * L a \ vX (x lP> q>, q) \ 

— Qstr< = >, (34) 

[ 2 % /2m c (fc2 - 2x 2P • q - 2q ■ k ± + ie) J 

where q' — q — x 2P — k 2 j_ is the momentum of outgoing quark after hard scattering, L a ^ i/X {xi P ,q l ,q) has 
been defined in Eq.(|l6|) and shown in Fig. 3. In Eq.(|34|) we have taken k 2 j_ = ~k 3 ±^ = k± by momentum 
conservation. In case of the antiquark participating the interaction, as shown in Fig. 5(b), we have 

L 2V2m c [(-q + x 2P + k 2± ) 2 -raj. + ie] pl J 

((~2p-q + ^)(4-mMQ,S z )(4 + m c))T^L a ^ x (x lP , q ,q')\ 

= Qs± T \ 7= 7, t ■ (00 J 

[ 2V2m c {k 2 x - 2x 2 p ■ q - 2q ■ k± + ie) j 

The total contribution of Fig. 5 is 



~9s 1 



2V2rn c , 2 + ^LZ^ ~ " 2 ' (36) 



2pq 



where 



LZl 2 Mp,k 2 ,k 3 ) = Tr 



m c MQ, + m c )(l - ^K)T^ il/X (x lP , q', q) 

z P ■ q 



-Ll\ vX (x lP , q, q')T a *(l ~ - m c )i{Q, S z )(4 + m c 



(37) 



In above Eq. ( |36| ) we have explicitly extracted the pole contribution from the quark or antiquark propagator 
in amplitude L^"\ (a) (xip, k 2 , h). 

Similarly, we can derive the corresponding result for amplitude -R^^ 4 4 from the diagrams shown in Fig. 6. 
Fig. 6(a) represents the case of soft gluon attaching to heavy quark line, which gives 

R asai(a) (x,v ko k.) - a Tr\ {i + ^ + k± + mc)M + TOc ^ Q ' Z m ^ T a*p** ( X1V a " ) 

((2p-q-jt ± M4 + rnMQ,S z )(4-rn c )T^Rll (q'', q )\ 

= 9sJ r < = - > , (38) 

( 2V2m c (kl + 2x 3P ■ q - 2q ■ k± - ie) J 

where q" = q + x 3P + k 3 ±. Fig. 6(b) shows the case of soft gluon attaching to heavy antiquark, it gives 



RlT; ib \x lP ,k 2 ,k 3 ) = q s Tr ( {i + m ^ SM - m ^~* ~ x ** - ^ ± ^ r*B%, u (x lP , q, q") 
A ^ 4 V M 3; ./ 1 2 ^2m c [{-q-x 3P -k 3i _y -ml-ie) Xl/ » iK ' 
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Figure 6: Graphical representation for the right amplitude -R^^ 4 4 in Eq. (p2[) for cc interacting with soft 
gluon, (a) shows the quark picks up the gluon, (b) shows the antiquark picks up the gluon 



g s Tr 



: + m c )j(Q, S z )(j + m c )(-2p ■ q + j>^)T- R a x l^{x lP , q, q") 
2\f2m c {k\ + 2xzV ■ q — 2q ■ k± — ie) 



(39) 



The total result of the amplitude i?^"* is 



9s 1 



? a 3 a 4 (f>) , 



2^/2: 



m, 



^3 



2q-k±-k 2 l 



(40) 



where 



2p-q 



){i + m c )j(Q,S z ){i-m c ) 



2p-q> 



: + m c )/(Q, S z )(4 - m c )(l - £^)T a 3 J R- t ( XlP , g , g ") 



(41) 



Now we can perform the integration through = 1, 2, 3) in function F{z^ , z^). From Eq. (pq) and 
(fi0|), we see that both variables xi and £3 have poles in the upper half complex plane. So after fixing x\ 
from 5 function, we can carry out the other integrals J dx^dx-i by using contour integrations in X2 and X3, 
closing at infinity, and circling the poles, we get 



F( z l > z 2 > z 3 ) — 



16-7T s + t — m| 8ml dk f [dk± p 



dxidx 2 dx^ XlP+ ^ +^P + z~+^x^+zl 
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1 r/ n-2Q-fc ± -fcj 



to + ^^S^ - fe) (-3 - " fe) « + *-»<;. 

x L^ A (a ;i p ) fc 2 ,fe3)®^(^iP»fc2,fe3)< I/14 5 0ia4 5 a2a3 } i n 
dQidy g 2 M(^ c ) d 2 



16-7T s + t — 8m 2 . dk^_dk± p 



J dx 2 dx 3 e ixiP+ Zl + ix zP + ( z 2 ~ z i )+ix 3 p + z~ 



(^2 + 2p . g - - IE) (S3 - l£ ) 



^ L££ (HP, fc 2 , k 3 ) ® (up, fc 2 , fc 3 ) ) 

J Oj_ 



ggjtfy g| ■M(V'c) O 2 ( e i x - L p+z-+ix 2 p+(z--x-)+ix 3 p+z- 

167r s + t — m% 8m 2 dk^dk^p X 
x (i27r(9(z 2 ~ - z 1 "))(i27r6»(z 3 "))^ lA14 ( 5 QlQ4 J a2a;i 

x Ll\ll(x lP ,k 2 ,k 3 )^RZ a ^(x lP ,k 2 ,k 3 ) _ a „. k± _ k , ± ) , (42) 

xi- X3- 2p - q ) k ± =0 ± 

where 

u-2Q-k ± -kl 

X\ = * — — . (43) 

s + t-ml y ' 

By substituting Eq.(E3) into Eq. (0o|) , we can obtain the lowest order real double scattering cross section 



dQ\dy 2s x 8 3 16tt(s + t - ml) 8m 2 J """ 2 2tt 



dz 1 dz 2 dz 3 —6(z z - z x )0{z 3 ) 



< p A \Ai{z^)G h + ^z 2 )G+^{z 3 )Al{0)\p A > d^5 a ^5™° 
d 2 



dk p L dk± p 



( e ix 1P +z 1 +ix 2P +(z 2 - Xl )+ix 3P +z 3 Ll\l\(x lP , k 2 ,k 3 ) 



2g-fcj_ — fc 

a^2 — —X3= 



(44) 



2p 9 J fc_i_=0 



One important step in getting the final result is taking the derivative with respect to k± as define in Eq. ( J27 
From Eq. (|37|) and (fill), it is obviously that when k± — 0, x 2 — x 3 — 0, and x% = x, therefore q' = q" = q, 
we found 

L^Vx(x lP , k 2 , k 3 )\ k±=0 = Rf v l\{x lP , k 2> k 3 )\ k±=0 = 0. (45) 

This indicates that there are cancellations in both the left and right amplitudes LV^\ and -Rjw* °f 
eq. (E3) . These are the cancellations between the amplitudes of gluon interacting with heavy quark shown in 
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Fig. 5(a), 6(a) and corresponding ones of antiquark cases in Fig. 5(b), 6(b). These eliminate the twist-two and 
even some twist- four terms contained in Eq.(|2^). Nevertheless certain twist four terms will survive and we 
study them in the following. 

Eq. (|45|) denotes that the derivatives on the exponential e M iP +z i +^ 2 p + {z 2 -z 1 )+ix 3 p+z 3 j nQ ^ con t r ibute, 
and that we can therefore set e lXlP+z i + tx ^p + ( z 2 ~ z i )+*x3p + z 3 _ e ixp z 1 _ ]\j ow we defined four gluon matrix 
element 

x <p A \A a± (z^)G b +± (z^)G b + x (z 3 )A a ± (Q)\p A >. (46) 
Substituting Eq.(ff6|) into Eq.(Q), we obtain the contribution from real double scattering as 

d<7 iA ] = _ 9 2 s 1 ^(x)^^) , 1W 3 

dQ\dy 2sx 8 3 8(s + t-ml) x 8m 2 c x 



& 



2 



{L^ll(x 1 p,k 2 ,k 3 )®RZl\(xiP,k 2 ,k 3 ) 2q . k± . h , | . (47) 

^ M = -I3 = 7: - 



2p g J k± =0 



Now we begin to consider the contribution from the virtual diagrams shown in Fig. 4(b) and Fig. 4(c). 
The graphical represents of the amplitudes corresponding to the left and right parts of diagram in Fig. 4(b) 
are given in Fig. 7 and Fig. 8 respectively, and the corresponding represents for Fig. 4(c) can also be given 
similarly. Following the same methods as that we use to derive Eq.(^?j), we obtain the differential cross 
section coming from virtual double scattering as 

d«T ,? i 



jA_ = _ _9t 1 M g (x) M^Pcl dtlll , igaiai§a2a3 

""'s + t — rrA) x 8ml ± 



dQ\dy 2s x 8 3 8(s + t~m^) x 8m 2 c 

9 



x 



2q-k ± -k z L 2q-k ± +k* 

X 2 ? n . n '^3 Onn 



+ ^vx^vMM) ® RlT l !" i {xiP,k 2 ,h) 2q . k± -ki 2q .k ±+ ki > , (48) 



2P '9 ■> k I =0 



where 



l a £T{xxPMM = Tr^T^{l-M-){i-m c )i{Q;S z ){i + m c ) 



( 1 + ^- q HU(^lP^^)\, (49) 
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(d) (c) 



fix, a\ /i2,a 2 /U3,a 3 




(a) (b) 

Figure 7: Graphical representation of the amplitude corresponding to the left part of diagram in Fig. 4(b); 
a) two soft gluons attach to quark line, b) two soft gluons attach to antiquark line, c) and d) one gluon 
attaches to quark line and another attaches to antiquark line 
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04, ^4 




Figure 8: Graphical representation of the amplitude corresponding to the right part of diagram in Fig. 4(b) 



2p-q' 



: + m c )^(Q,S z )(i- m c ) 



(1-^)^(^,9X^2)^. 



TV 



m c )i{Q, S z )(4 + m c )L a i ^ vX {xp, q, q) 



^At M4 ( x iP' fc 2,A:3) = Tr (4 + m c )^(Q,S z )(4-m c )R'^ IJ/i (xp,q,q) 



(50) 

(51) 
(52) 



where 



xi 

Q2 



q-p 

q - x 2 p - kx , 
q - x 3 p + k± . 



Now we can work out the derivative, the calculation is tedious but straightforward. The final result for 
lowest order double scattering cross section is 



dQ\dy 



da 



dQ\dy dQ\dy 
16TT 3 aa 3 s e 2 Q \R(0)\ 2 M g/A (x) 



f{s,i ml) 



(S-m^(i-m^(S + ty : 



(53) 
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where the function / is 



/(s, t, ml) = {-72Q 2 ± m$s 2 + 16Q 2 L mfs 2 {82s + 73*) + 8Q 2 ± m$ s 2 {-179s 2 - 550 si - 88* 2 ) 

+ 8Q 2 L mlli>s 2 {-89s 3 + 518s 2 * - 135s* 2 - 158?) 

+ 8Q\m%s 2 {Z87s i + 290s 3 * + 1327s 2 * 2 + 1204s* 3 + 337* 4 ) 

+ 8Q 2 l to^,s 2 (-402s 5 - 870s 4 * - 2263s 3 * 2 - 2560s 2 ? - 1811s* 4 - 516* 5 ) 

+ 8Q 2 ± mp 2 (128s 6 + 498s 5 * + 1529s 4 * 2 + 2222s 3 * 3 + 2182s 2 * 4 + 1160s* 5 + 306* 6 ) 

+ 256Q 2 L m^,s 3 t(-s 5 - 3s 4 * - 5s 3 * 2 - 5s 2 * 3 - 3s* 4 - * 5 ) 

+ 9mJ, 8 (s 2 + 2s* + 1 2 ) + 4m 1 / 5 (-53s 3 - 23s 2 * + 17s* 2 - 13* 3 ) 

+ mj, 4 (-181s 4 - 788s 3 * - 1552s 2 * 2 - 1208s* 3 - 263* 4 ) 

+ m.J 2 (1081s 5 + 3875s 4 * + 6582s 3 * 2 + 7030s 2 * 3 + 4537s* 4 + 1295* 5 ) 

+ m,J,°(-1296s 6 - 5283s 5 * - 10856s 4 * 2 - 14170s 3 * 3 - 12382s 2 * 4 - 7039s* 5 - 1958* 6 ) 

+ tt4(823s 7 + 2921s 6 * + 6357s 5 * 2 + 11051s 4 * 3 + 13065s 3 * 4 + 9643s 2 * 5 + 4659s* 6 + 1289* 7 ) 

+ m^(-224s 8 - 615s 7 * - 525s 6 * 2 - 1609s 5 * 3 - 4268s 4 * 4 - 4537s 3 * 5 - 2267s 2 * 6 - 843s* 7 - 320* 8 ) 

+ mis*(-36s 7 - 727s 6 * - 1561s 5 * 2 - 894s 4 * 3 + 290s 3 * 4 + 117s 2 * 5 - 389s* 6 - 192* 7 ) 



4 Numerical results and discussion 

For the purposes of numerical evaluation, we evaluate the ratio of transverse momentum double scattering 
differential cross section to corresponding single scattering contribution 



+ 3m|s 3 * 2 (44s 5 + 173s 4 * + 239s 3 * 2 + 116s 2 * 3 - 13s* 4 - 19* 5 )} /9. 



(54) 



In deriving (|53|) we have used the following relation for the long distance matrix element 



<0|Of( 3 5 1 )|0>=^|i?(0)| 2 . 



(55) 





(56) 
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We perform the integration through rapidity y in Eq.(|l9|) and Eq.(f33|) and the integration range are taken 
\y\ < 0.4. Before embarking on the numerical calculation, we have to know the non-perturbative parton 
distribution f g /A( x ) an d multi-parton correlation function M g (x). 

In Eq.(|ll]), the effective nuclear gluon distribution function f g /A{ x ) should have the same operator dcfin- 
ion of the normal parton distribution with free nucleon states replaced by the nuclear states. For a nucleus 
with atomic number A, we define 

fg/A{x) = Af g/N (x), (57) 

where / ff /jv(a;) is normal gluon distributions in a free nucleon. In this effective nuclear parton distribution, 
we have neglected the EMC effect and the difference between the normal parton distribution in a free neutron 
and in a free proton. We adopt f g /N{ x ) from 

The authors of JD| [l6| proposed following phenomenological expressions for the twist-four matrix elements 
that occur in soft rescattering: 

M i (x) = A 1 / 3 X 2 f i/A (x) ) (58) 

where i = q,q, and g. The fi/A( x ) ar e the effective twist-2 parton distributions in nucleus defined in Eq.(|57]). 
The constant A 2 has units of mass squared due to the dimension difference between twist-4 and twist-2 matrix 
elements. The A 2 should be determined by experimental measurement. It was estimated in fl3|| by using 
the measured nuclear enhancement of the momentum imbalance of two jets in photon-nucleus collisions, and 
was found to be order of 

A 2 re 0.05- O.lGeV 2 . 
In our following calculation we use A 2 = 0.05GeV 2 . 

From the definion of the correlation functions in Eq.ji^), the lack of oscillation factors for both z% and 
integrals can in principle give nuclear dependence proportional to A 2 I Z . The A 1 / 3 dependence is a result 

of the assumption that the positions of two field strengths (at z^ and z^ , respectively) are confined within 

one nucleon. 

In Fig. 9 we show the ratio of double scattering and single scattering versus Q±. We plot pW as a function 
of Q± for s = (20GeV) 2 and s — (40GeV) 2 respectively. The ratio is normalized by A 1 / 3 . We find that the 
overall contribution of double scattering give negative sign, while it is positive in the processes of Drell-Yan 
and DIS. From Fig. 9, it is obviously that the absolute value of the ratio increase with Q±, this indicates that 
in high Q± J/i[> is more suppressed by the double scattering effect than that in low Q±. We find that the 
two curves for s — (20GeV) 2 and s — (40GeT^) 2 are very close, this indicates that at this range the ratio is 
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2 4 6 8 10 12 

Qj_(GeV) 

Figure 9: The ratio of double scattering contribution and single scattering contribution normalized by A 1 / 3 
for s = 20GeV and s = AOGeV. 

almost independent of the invariant s. Although our quantity result is not large, about 17% at Q± = 8GeV 
for s = 40GeV^, but if we include the extra A 1 ! 3 factor, it can give observable effect in large A case. More 
over our result lies in reducing some theoretical uncertainties. 

In summary, we have applied LQS method, which was developed in inclusive DIS process, to study the 
multi-scattering effect on J/ip photoproduction. We consider the interaction between nuclear matter and 
cc pair produced from the hard collision of a photon and a parton in nucleus. This interaction undergoes 
through either quark or antiquark absorbing a soft gluon. We find the lowest order double scattering has not 
a great influence and only contributes a slight suppression comparing with the single scattering cross section. 
Our analysis can be extended to study the multi-scattering effects in J/ip production from hard-nucleus or 
nucleus-nucleus collisions, which will involve initial state strong interactions that are absent in the case of 
photoproduction. These two hadron production processes will be more complicated because for initial state 
we must consider two types of double scattering, soft-hard scattering and double hard scattering. 
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